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ABSTRACT 

The  air-bubbling  technique  utilized  by  the  Military  Sea 
Transportation  Service  to  prevent  ice  formation  in  North  Star 
Bugt,  Thule,  Greenland  is  discussed.  A  proposed  model  of 
Induced  water  circulation  is  presented.  Physical  processes 
impeding  ice  formation  and  growth  in  sea  water  are  described. 
Oceanographic  data  collected  in  conjunction  with  the  MSTS 
project  are  analyzed  and  presented  in  the  appendixes. 


FORE\roRD 


The  formation  of  sea  ice  in  northern  waters  often  hastens  the 
termination  of  shipping  at  Arctic  harbors  and  sites.  Retardation  of 
ice  growth  can  reduce  or  prevent  damage  to  marine  structures  such  as 
De  Long  pier  at  Thule,  Greenland.  A  thorough  understanding  of  physi¬ 
cal  effects  which  delay  ice  formation  and  slow  ice  growth  is  required. 
For  these  reasons,  the  experiment  performed  at  Thule  in  1959  is  of 
considerable  interest. 

This  report  is  a  study  of  the  effects  of  air-bubbling  on  the  physi¬ 
cal  properties  of  the  water  adjacent  to  De  Long  pier.  It  attempts  to 
formulate  a  working  hypothesis  for  explaining  the  mechanism  of  the  proc¬ 
esses  which  retarded  formation  and  growth  of  sea  ice. 

Conclusions  expressed  in  this  report  may  require  revisions  as 
additional  data  become  available.  All  additional  information  which 
might  amplii'y  or  modify  this  report  will  be  welcomed  by  the  I^dro- 
graphlc  Office. 
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INTRODUCTION 


Successful  application  of  an  air  bubbllnc  technique  for  prevent¬ 
ing  Ice  formation  during  the  Tall  of  195^^  enabled  the  Military  Sea 
Transportation  Sei'vlcc  to  c:d;end  the  shipplnf?  season  at  Thule,  Green¬ 
land,  The  normal  shipping  season  extends  from  early  July  to  early 
October.  Shipping  during  the  first  half  of  July  Is  usually  dependent 
on  Icebreaker  escort;  shipping  Is  ordinarily  terminated  prior  to 
Initial  ice  formation  In  autumn. 

Atlaptatlon  of  a  method  origlnaLly  developed  In  Scandinavian  coun¬ 
tries  for  prevention  of  fresh-water  freezing  permitted  maintenance  of 
an  Ice-free  area  (polynya)  adjacent  to  De  Long  pier  despite  normal 
ice  formation  In  the  surrounding  waters  of  North  Star  Bugt  (Bay). 
Safeguarded  against  becoming  frozen-ln  at  the  pier,  ships  of  the 
supply  convoy  remained  at  Thule  until  25  October  195Q  -  the  latest 
date  on  which  MSTS  had  ever  operated  In  such  a  northerly  location. 
Success  of  the  temporary  installation  prompted  the  Comnander  Military 
Sea  Transportation  Service  Atlantic  (COMSTSLANT)  to  formulate  plans 
for  the  establishment  of  a  permanently  installed  bubbling  system  at 
De  Long  pier. 

In  the  fall  of  1959#  the  U.  S,  Navy  Iftrdrographlc  Office  was 
requested  to  conduct  oceanographic  studies  concurrently  with  the 
operation  of  the  system  In  order  to  obtain  Information  on  physical 
processes  impeding  Ice  foi*matlon  and  growth  in  the  bay.  The  over¬ 
all  operation  was  dubbed  "Project  Polynya", 

NORTH  STAR  BUGT 

North  Star  Bugt,  approximately  3  square  miles  In  area,  recedes 
about  1-1/2  miles  northeastward  between  Astro  Pynt  and  Mount  Dundas 
on  the  southern  shore  of  Wolstenholme  Fjord.  The  entrance  of  the 
bay,  about  3/^  mile  wide,  is  narrowed  by  De  Long  pier  and  a  causeway 
which  extend  approximately  O.U  mile  west -northwestward  from  Astro 
Pynt*  The  pier,  1,000  feet  In  length  and  50  feet  In  width,  adjoins 
the  causeway  and  Is  paroUel  to  It,  Inside  the  entrance,  the  width 
of  the  bay  Increases  to  about  1  mile. 

POLYNYA  INSTALLATION  1959 

Under  the  supervision  of  MSTSLANT,  Canadian  Uhderwater  Demolition 
Unit  BRAVO  began  Installation  of  the  bubbling  system  In  North  Star  Bugt 
early  In  September  1959*  Briefly,  the  system  consists  of  perforated 
submarine  polyethylene  pipes  which  serve  to  conduct  compressed  air  to 
the  bottom  of  the  bay  and  distribute  It  in  the  form  of  bubbles  over  a 
wide  area  (see  Figure  1).  An  auxiliary  alcohol- Injection  system  Is 
utilized  to  prevent  or  eliminate  ice  formation  within  the  pipes. 

Installation  was  completed  on  27  September,  Equipment  was  Inter¬ 
mittently  tested  until  appreciable  ice  began  to  form  In  the  bay  on 
7  October,  at  which  time  full  operation  was  begun. 
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ICE  FORMATION  AND  GROWTH  IN  1959 


Ice  was  first  observed  on  l6  September  along  the  eastern  fringe 
of  North  Star  Bugt  at  a  point  where  fresh  water  flows  from  the  Pltuflk 
River*  Hy  24  September  grease  ice  formed  In  an  area  northeast  of  the 
pier,  where  shallow  water  was  observed  to  cool  to  the  freeelng  point 
during  low  tide.  This  grease  ice  drifted  westward  with  the  next  ebb 
tide.  Ry  7  October  a  considerable  amount  of  pancake  ice  had  formed  in 
the  area  northeast  of  the  pier. 

During  the  morning  of  8  October  a  foehn  wind  arose  with  gusts  to 
approximately  50  knots.  Ihe  air  temperature  Increased  from  17°  to 
34°  F*  On  the  following  morning,  the  bay  was  completely  clear  of  Ice 
and  the  surface  water  temperature  had  Increased  from -1*22”  to  -0*94°  C* 

During  the  morning  of  10  October  grease  and  pancake  ice  of  small 
■diameter  formed  over  most  of  the  area  northeast  of  the  pier.  On  11 
October  a  considerable  amount  of  large  pancake  ice  was  observed  drift¬ 
ing  Into  the  region  from  the  south.  Hie  areas  to  the  south,  west,  and 
northwest  of  the  pier  attained  a  coverage  of  approximately  nine-tenths) 
no  major  Ice  formation  was  observed  near  the  pier*  Ry  15  October  most 
of  the  bay  was  covered  with  young  floes  averaging  approximately  10  feet 
in  diameter.  Ten-tenths  concentration  of  young  Ice  accumulated  every¬ 
where  by  20  October' except  In  the  Ice-free  area  adjacent  to  the  pier* 

DATA  COLLECTION 

Serial  temperature,  salinity,  and  current  data  were  obtained  at 
the  locations  shown  In  Figure  2*  Temperature  and  salinity  observa¬ 
tions  were  taken  between  9  September  and  21  October* 

Station  1  was  occv^led  dally  using  a  small  hand  winch  mounted  on 
the  pier  or  by  use  of  an  oceanographic  winch  aboard  the  VESIVIND.  For 
detailed  observations,  the  polynya  was  divided  Into  Stations  lA,  ID, 
and  1C  as  shown  In  Figure  3*  Station  ID  was  designated  at  a  point  ly¬ 
ing  approximately  15  feet  outside  the  polynya  area  and  near  the  east¬ 
ern  end  of  the  pier*  Stations  2,  3>  end  4  were  occupied  twice  weekly 
using  either  a  Greenland  cruiser  or  an  LCVP.  Station  5  occupied 
only  on  15  September  and  Station  ID  only  on  13  October*  Stations  6, 

7,  8,  and  9  eere  occupied  weekly  by  the  WESIVIND. 

A  200-foot  bathythermograph  was  used  to  measure  temperature  at 
all  oceanographic  stations*  RT  drops  were  made  daily  within  the 
polynya  area*  Original  plans  had  Included  dally  observations  at  Sta¬ 
tions  1,  lA,  and  16;  however,  presence  of  shipping  often  precluded 
data-collection  at  all  3  locations*  A  bucket  thermometer  was  used  In 
conjunction  with  each  BT  drop  to  obtain  siU'face  water  temperature* 
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Meteorological  observations  obtained  dally  on  the  pier  Include 
wind  speed  and  direction,  wet  and  dry  bulb  air  temperature,  cloud 
cover  and  type,  visibility,  sea  state,  and  state  of  weather.  In 
addition,  pertinent  data  were  extracted  from  the  weather  log  main¬ 
tained  at  Thule  Air  Base  for  analysis  (Appendix  V).  Data  on  tidal 
currents  were  obtained  by  suspending  3  Roberts  radio  current  meters 
at  depths  of  6,  26,  and  5^  feet  from  an  anchored  radio  buoy.  Water 
depth  at  each  location  was  12  fathoms.  Signals  transmitted  from 
the  buoy  were  recorded  at  30-nilnute  Intervals  by  a  monitoring  sta¬ 
tion  aboard  the  WESIWIND.  Current  meter  Station  10  was  occupied  from 
I83OZ,  25  September  to  2100Z,  26  September;  Station  11  was  occupied 
from  1200Z,  6  October  to  1930Z,  13  October  1959*  Current  data  were 
not  tabulated,  because  the  recorded  results  either  approximated  the 
threshold  value  of  the  current  meter  (0.2  knot)  or  were  unreadable. 

A  supplementary  survey  conducted  during  April  i960  by  Hydro¬ 
graphic  Office  ice  observers  yielded  late-winter  temperature  data 
at  Stations  1  and  2  (Appendix  III). 

DATA  ANALYSIS 

Oceanographic  conditions  at  each  station  were  examined  for 
factors  contributing  to  the  formation  and  growth  of  ice.  Surface 
temperatures  Indicating  heat  lose  at  the  sea  svirface  and  pljysical 
properties  showing  the  distribution  of  heat  loss  throughout  the 
water  coliunn  were  studied.  Data  obtained  outside  the  polynya 
were  compared  to  those  obtained  at  Station  1  in  order  to  deter¬ 
mine  the  effect  of  the  bubbling  system  on  the  oceanographic  struc¬ 
ture. 


Reversal  of  the  heat  budget  had  occurred  prior  to  Inception 
of  the  oceanographic  survey.  Except  for  Interruption  by  the  foehn 
on  8  October,  progressive  cooling  was  observed  at  sdl  depths.  The 
temperature  rise  shown  by  the  tabulated  data  for  9,  10,  and  11  Octo¬ 
ber  was  observed  at  Station  1.  Upon  resumption  of  the  cooling  pro¬ 
cess,  surface  temperature  outside  the  bubbled  area  decreased  rapidly. 
The  freezlxig  point  was  attained  on  13  October. 

A  study  of  the  salinity  structure  Indicates  spatial  and  tempo¬ 
ral  fluctuations  of  surface  values  and  depth  of  the  Isohallne  layer. 

A  plot  of  the  surface  salinity  values  at  Station  1  Is  presented  In 
Figure  4.  Hie  portion  of  the  plot  constructed  from  values  for  early 
September  Indicates  that  a  certain  periodicity  may  exist.  The  pro¬ 
nounced  Increase  during  the  latter  part  of  September  Is  attributed 
to  cessation  of  runoff. 

Data  obtained  at  Station  3  on  18  September  and  12  October  are 
plotted  in  Flgiures  ^  and  6.  The  surface  water  temperature  on  18 
September  was  0.64°  C;  temperature  maxlmvon  of  0.72°  C  occurred  at 
9  and  20  meters.  The  surface  salinity  was  31 *23  °/oo.  Convection 
extended  to  a  depth  of  only  a  few  meters. 
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Sy  12  CXstober  the  siirface  vater  had  cooled  to  -1*32°  C;  the  vartn- 
est  water  was  at  the  hotton*  Ihe  temperature  maximuin  of  -0.62°  C  was 
observed  at  a  depth  of  33  meters,  mie  surface  salinity  had  Increased 
to  32.34  ^/oot  and  convective  mixing  had  produced  an  Isohallne  layer 
in  the  upper  1$  meters.  Die  calculated  freezing  point  of  the  surface 
is  -1.76°  C. 

Hie  density  gradient  below  the  13-metcr  level,  although  weak,  has 
significant  relevance  to  the  bubbling  system.  O^oretlcal  ice-poten- 
tlal  calculations  using  the  data  of  the  deeper  stations  show  that, 
prior  to  initial  ice  formation  »  thermohallne  convection  takes  place 
to  a  depth  of  approximately  15  meters.  Consequently,  the  water  below 
this  level  temporarily  serves  as  a  source  of  sensible  heat.  However, 
as  ice  forms,  the  salinity  of  the  upper  layer  increases,  resulting  in 
greater  density  and  an  increase  in  depth  of  convective  mixing.  As  the 
density  gradient  weakens  and  eventually  disappears,  cooling  to  the 
freezing  point  will  occur  throughout  the  water  column.  At  the  known 
average  rate  of  heat  loss  from  the  sea  surface  in  the  latitude  of  ^niule, 
the  entire  svqpply  of  warm  water  will  be  eliminated  within  two  weeks 
after  initial  formation  of  ice. 

The  heat  content  of  air  Issuing  from  the  compressors  is  consider¬ 
able.  An  appreciable  amount  of  beat  is  possibly  introduced  into  the 
bottom  water  adjacent  to  the  pier  when  compressed  air  cools  in  the 
polyethylene  pipes;  however,  the  data  do  not  indicate  a  tenqoeirature 
differential  attributable  to  this  source  within  the  bubbler  field.  A 
layer  of  dirt  covering  four  steel  feed  pipes  provides  insulation;  how¬ 
ever,  heat  loss  through  the  rubber  feed  hoses  is  great.  Between  the 
point  where  these  hoses  connect  to  the  steel  pipes  and  the  point  where 
they  enter  the  water,  melting  of  snow  within  a  radius  of  approximately 
2  feet  was  observed.  Ice  formation  due  to  moisture  condensate  in  the 
underwater  sections  of  the  feed  pipes  was  removed  by  alcohol  injection. 

The  effectiveness  of  the  bubbler  system,  when  warmer  bottom  water 
is  available,  is  manifested  by  the  temperature  data  in  Ibble  I. 


TABLE  I 

13  October  1959 


STATION  1  STATION  ID 


Depth 

Temp. 

Depth 

Temp. 

(meters) 

(°c) 

(meters) 

(°c5 

0.0 

-1.09 

0 

-1.76 

3.5 

-1.08 

3 

-1.60 

6.5 

-1.06 

6 

-1.25 

9.5 

-1.06 

9 

-1.24 
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The  surface  temperature  at  Station  ID,  located  Imnedlately  outside 
the  bubbler  field,  shows  that  the  surface  water  had  cooled  to  the 
freezing  point.  Despite  ice  formation  around  the  perimeter  of  the  agita¬ 
ted  area,  the  data  obtained  at  Station  1  show  the  surface  temperature  to 
be  0.67°  C  above  the  freezing  point.  The  data  for  Station  ID  are  assvnned 
to  be  indicative  of  temperature  data  that  would  have  been  observed  at 
Station  1  had  the  bubbler  system  not  been  in  operation.  Ttoe  data  of  13 
October  plus  the  profiles  for  Stations  2  and  3  on  12  October  Indicate 
that  water  from  depths  greater  than  15  meters  is  circvLLated  into  the 
agitated  water  columns  adjacent  to  the  pier. 

Proof  that  the  bubble  system  acts  as  a  huge  pump  capable  of  perform¬ 
ing  work  on  the  surroxinding  subsurface  water  is  provided  by  comparison 
of  data  presented  in  Tables  II  and  III. 


TABLE  II 
11  October  1959 

STATION  1 


Depth 

(meters) 

Temp. 

(°  c) 

Salinity 

(%o) 

0.0 

-0.90 

32.^2 

26.08 

3.5 

-0.91 

32.il2 

26.08 

6.5 

-0.91 

32.42 

26.08 

9.5 

-0.89 

32.43 

26.09 

TABLE  III 
12  October  1959 

STATION  2 


Depth 

(meters) 

of 

Salinity 

(o/oo) 

0 

-1.42 

32.32 

26.02 

5 

-1.40 

32.32 

26.02 

15 

-1.42 

32.32 

26.02 

22 

-0.69 

32.52 

26.16 

24 

-0.78 

32.52 

26.16 
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The  density  of  the  agitated  water  column  In  Table  II  is  greater 
than  the  density  to  at  least  1^  meters  In  Table  III;  therefore,  work 
was  performed  by  the  system  In  raising  water  through  a  -vertical  dis¬ 
tance  In  excess  of  13  meters.  Ccxnparison  of  salinity  and  density 
data  of  Table  II  and  the  plotted  curves  of  Figure  6  reveals  that 
water  similar  to  the  entire  water  column  at  Station  1  is  found  at  20 
meters  at  Station  3>  indicating  that  the  water  was  raised  at  least 
20  meters. 

The  eventvial  cooling  of  the  entire  water  column  to  the  freezing 
point  indicates  that  vertical  transport  of  sensible  heat  from  depth 
was  not  a  factor  in  the  maintenance  of  the  artificially  created 
pdlynya,  except  during  the  initial  stage  of  the  ice  formation.  Con¬ 
sequently,  an  understanding  of  the  physical  process  Involved  must  be 
sought  along  other  lines. 

Elementary  Ice  particles  are  probably  disk-shaped  and  devoid  of 
crystalline  form.  Ordinarily  they  flocculate  and  grow  into  true 
cxystals.  The  turbulent  energy  of  the  induced  currents  may  destroy 
the  crystals  before  they  enlarge  or  may  effectively  prevent  crystal¬ 
line  growth  about  Ice  nuclei.  Ice  particles  at  the  surface  of  the 
bubbled  area  are  rapidly  swept  from  regions  of  divergence  into 
regions  of  convergence  where,  by  means  of  descending  currents,  they 
are  transported  beneath  the  surface  to  be  eventually  dispersed  from 
the  polynya  area. 

lydrographic  Office  ice  observers,  stationed  at  Thule  Air  Base 
throughout  the  winter  of  1959-60,  noted  that  the  polynya  gradually 
narrowed;  by  the  end  of  December  width  ranged  from  12  feet  at  the 
eastern  end  to  50  feet  at  the  western  end,  where  an  auxiliary  air 
hose  was  used  to  augment  the  bubbling  activity  by  inducing  more  vig¬ 
orous  currents.  Dimensions  of  the  ice-free  area  gradually  increased 
during  spring  as  the  air  temperature  rose  to  approximately  0°  F. 

A  plot  of  sea  ice  tensile  strength  versus  temperature  (Assur, 
1956)  shows  a  marked  increase  of  strength  as  the  temperature  of  the 
ice  drops  below  -9.2°  F.  At  this  temperature  sodium  chloride  is 
precipitated  from  the  brine  pockets  in  the  ice.  During  periods  of 
extremely  low  air  temperature  in  winter,  the  weakest  point  of  the 
ice  should  be  at  its  undersurface  where  the  temperature  approaches 
that  of  the  water. 

Measurements  made  during  April  i960  show  that  ice  thickness 
directly  above  one  of  the  polyethylene  pipes  averaged  approximately 
10  inches  while  thicknesses  ranged  between  4l  and  44  inches  at 
locations  60,  200,  and  379  yards  north- northeast  of  the  pier.  Abra¬ 
sive  action  of  induced  currents  apparently  inhibited  ice  growth  in 
zones  of  most  vigorous  flow. 
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The  erosive  capability  of  water  currents  is  manifested  by  recent 
experiments  in  the  Antarctic,  Specially  shaped  propellers  driven  by 
small  motors  were  suspended  through  holes  in  the  ice  of  McMurdo  Sound. 
The  propellers  created  vigorous  currents  which  eroded  the  ice  from 
below.  A  10-horsepower  device  reportedly  required  I83  hours  to  open 
an  area  30  by  85  feet  in  8-foot-thick  ice.  An  additional  swath  of  ice 
200  feet  long  was  eroded  to  a  thickness  of  I8  inches;  soon  afterward, 
it  fell  through. 

Analysis  of  data  obtained  with  the  Roberts  current  meters  revealed 
no  permanent  current.  Mass  transport  of  water  in  the  area  was  attri¬ 
buted  to  tidal  action.  Peak  tidal  current  speed  was  approximately  0.2 
knots  (based  on  threshold  value  of  the  Instrument). 

CONCLUSIONS 

The  bubbling  system  operates  as  a  huge  pump  capable  of  performing 
work  on  contiguous  subsurface  water.  The  rising  streams  of  bubbles 
initiate  a  system  of  circulatory  cells  which  extend  from  the  bubbled 
region  into  adjacent  water.  Water  from  depths  exceeding  15  meters  is 
circulated  into  the  agitated  columns  adjacent  to  the  pier  and  brought 
to  the  surface. 

At  the  time  of  initial  ice  formation  in  1959;  convective  mixing 
had  occurred  throughout  the  upper  15  meters  of  North  Star  Bugt.  The 
density  gradient  below  the  15~roeter  level  gradually  weakened  with  ice 
growth,  and  the  entire  water  column  cooled  to  the  freezing  point. 

After  elimination  of  the  warm  water  supply,  maintenance  of  an  open 
water  area  adjacent  to  the  pier  was  attributed  to  the  ice-dispersive 
and  erosive  activity  of  the  induced  currents  coupled  with  the  possi¬ 
bility  that  the  turbulent  energy  also  sufficed  to  prevent  crystalline 
growth  about  ice  nuclei.  Efficiency  of  the  system  varied  directly 
with  turbulence. 

Considerable  nurrowini;  of  the  polynya  Ly  mid-winter  was  attrib¬ 
uted  to  marked  increase  of  terisile  streni^th  with  consequent  Increased 
resistance  to  erosive  action  of  the  induced  currents  as  the  tempera¬ 
ture  of  the  ice  dropped  below  -9.2°  F.  Vertical  growth  of  the  newly 
formed  ice  cover  within  the  buooled  area  was  inhibited  by  this  ero¬ 
sive  action  because  the  undersurface  of  the  ice  is  weakest  wnen  its 
temperature  equals  that  of  the  water.  Increase  in  the  dimensions  of 
the  ice-free  area  was  observed  to  concur  with  an  increase  of  air 
temperature  to  approximately  0°  F  in  early  spring.  This  increase  was 
attributed  to  marked  decrease  of  tensile  strength  with  consequent 
decreased  resistance  to  erosion  as  the  temperature  of  the  ice  rose 
above  -9.2°  F. 
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Uhlque  properties  of  fresh  water  make  the  bubbling  system  highly 
suitable  for  lakes  and  to  a  somewhat  lesser  extent  for  brackish  estu¬ 
aries.  The  system  Is  less  effective  In  salt  water,  because  maximum 
density  of  water  with  salinity  in  excess  of  2U,7  o/oo  Is  attained  at 
the  freezing  point.  However,  factors  other  than  the  upward  transport 
of  warm  water,  as  previously  discussed,  also  contribute  to  the  maln- 
enance  of  an  Ice -free  area. 

In  regions  where  upward  circulation  of  sensible  heat  Is  not  a 
factor,  maintenance  of  an  ice-free  area  is  predominantly  dependent 
upon  speed  and  Intensity  of  the  Induced  currents. 
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FIGURE  5  TEMPERATURE ,  SALINITY,  AND  DENSITY 
PROFILES  AT  STATION  3,  18  SEPTEMBER 
1959 


FIGURE  6  TEMPERATURE,  SALINITY,  AND  DENSITY 
PROFILES  AT  STATION  3,  12  OCTOBER 
1959 


BIBLIOGRAPHY 


ASSUR,  A.  Coii5)08itlon  of  sea  ice  and  Its  tensile  strength,  National 
Academy  of  Sciences  -  National  Research  Council  Publication 
271  P»  December  195^* 

DOTSON,  W.  A.  and  McINROY,  G.  I.  Project  reports;  Operation  Polynya, 
1958-1960,  [n.  p.2  Unpublished. 

LAMB,  H.  Hychrodynamlcs.  New  York:  Dover  Publications.  730  p.  19^5. 

LEE,  0.  S.  and  SIMPSON,  L.  S.  A  practical  method  of  predicting  sea 
Ice  formation  and  growth,  U.  S.  Navy  Iferdrographlc  Office  Technical 
Report  No.  U,  27  p.  September  195^. 

McGOUGH,  R.  J.  Local  environmental  factors  affecting  ice  formation 
in  North  Star  Bugt,  Greenland,  U.  S.  Navy  I^drographlc  Office 
^technical  Report  No.  23,  **0  p.  January  195o. 

U.  S.  AIR  FORCE  dally  weather  records;  Thule,  Greenland,  1959>  19^0. 
Unpublished. 


13 


APPENDIX  I 


MODEL  OF  THE  POLYNYA 
CIRCULATORY  SYSTEM 


APPErJDIX  I 


MODEL  OF  THf:  POLYWA  CIRCULATORY  SYSTEM 


Flytire  7  llluctrates  streamlines  ncnerated  by  the  motion  of  a 
solid  sphere  in  an  infinite  mass  of  frictionlecs  fluid.  If  we  take  the 
origin  at  the  center  of  the  sphere  and  the  x  axis  In  the  direction  of 
motion,  the  normal  velocity  (Vp)  at  tne  suivrnce  of  the  sphere  is  U  cos  0, 
where  U  is  tlie  velocity  of  the  center. 


Lamb  shows  that  the  stream  function  duo  to  the  sphere  is 

sin^e 


whore  a  is  the  radius  of  the  sphere,  r  is  the  radius  vector  from  the 
center  to  points  on  or  exterior  to  the  sphere  (r*a),  and  d  is  the  ancle 
between  the  radius  vector  and  the  x  axis.  At  any  civen  instant  the 
trajectories  of  the  fluid  particles  are  tancent  to  the  streamlines. 

The  total  flux  through  a  curved  surface  S  is 
Arbitrarily  making  this  value  equal  to  -2Tri^,  we  have 


-2ir 


In  the  case  where  S  is  the  surface  of  the  above  sphere  (rsa) 
tion  of  2'nryd3  for  dS  yields 


substitu- 


-27rV^=j^Vn27rydS, 


where  ds,  as  shown  in  [‘Ic.urc.  i ,  is  an  infinitesimal  length  of  arc  sub¬ 
tended  by  an  inf initcsmal  ariijuc,  d  9  ,  on  the  sur  'aco  S.  Substitution 
of  U  cos  9  ,  a  sin  9  ,  and  a  dfi  for  Vp,  y,  and  ds,  respectively,  and 
Intcgratinf]  between  the  lii.iits  0  and  9  yields 


Therefore, 


^z-lUa^sin^  9. 


Iy.»mb  shows  that  the  stream  I’unction  from  an  n  pole  is  given  by 

a""'  a 

OK 

Since  the  sphere  acts  as  a  dipole, 

dx  r  • 

From  the  boundary  value  rsa, 

^8in*  fl=--5Ua^  sin^  9. 

0  z 
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Therefore, 


K  =  --^Ua® 

2 

and  for  the  general  cacej  i.  e.,  r>0 

A  model  of  the  polynya  circulatory  system  can  be  formulated  from 
the  idealized  case  by  adaptation  of  the  principles  to  the  bubbling  sys¬ 
tem,  Considering  the  motion  of  each  ascending  bubble  to  be  directed 
along  the  positive-downward  Z  axis,  there  will  be  a  streamline  coinci¬ 
dent  with  the  Z  axis  and  a  vertical  flow  of  water  particles.  Ascending 
motion,  represented  by  a  negative  vertical  velocity,  creates  divergence 
at  the  surface.  Approximately  midway  between  bubble  streams  is  a 
region  of  convergence  with  consequent  descending  motion,  clearly  dis¬ 
cernible  in  Figure  9. 

Surface  water  beyond  the  pipe  furthest  from  the  pier  flows  out¬ 
ward  to  a  distance  determined  by  the  horizontal  momentum  of  the  water 
particles. 

The  data  show  greater  density  in  water  brought  to  the  surface  by 
the  bubble  activity  during  the  pre-freezeup  and  initial  freezeup  peri¬ 
ods.  Consequently,  as  the  higher  density  surface  water  flowing  outward 
from  the  divergence  zone  above  pipe  #4  suffers  a  gradual  decrease  in 
the  horizontal  component  of  the  velocity  vector,  the  vertical  component 
increases.  From  the  point  where  the  horizontal  component  becomes  zero, 
descending  motion  extends  to  depths  where  divergence  directs  a  horizon¬ 
tal  component  toward  the  pier. 

The  proposed  model  of  the  polynya  circulatory  system  is  presented 
in  Figure  10.  This  cross-sectional  view  shows  the  eastern  ends  of  the 
polyethylene  pipes;  arrows  Indicate  principal  paths  of  the  water  parti¬ 
cles. 


The  author  is  indebted  to  Dr.  Lloyd  Simpson  of  the  Iferdrographlc 
Office  for  advice  and  assistance  in  application  of  hydrodynamic  prin¬ 
ciples  in  development  of  this  idealized  model  of  the  bubbling  system. 
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FIGURE  9  FLOW  PATTERN  AT  THE  POLYNYA  SURFACE 


FIGURE  10  VERTICAL  CROSS  SECTION  SHOWING  QRCULATION  PRODUCED  BY  BUBBLING  SYSTEM 
IN  NORTH  STAR  BUGT 


APPENDIX  II 


OCEANOGRAPHIC  PROGRAM  -  I960 


APPENDIX  II 


OCEANOGRAPHIC  PROGRAM  -  I96O 

A  program  of  oceanographic  data  collection  similar  to  that  of 
1959  was  conducted  during  the  fall  of  i960.  Observations  of  the 
formation  and  growth  of  Ice  on  North  Star  Bufit  were  initiated  on  6 
October.  On  this  date  grease  and  pancake  ice  were  observed  in  the 
shallow  water  area  northeast  of  De  Long  pier.  By  30  October  most 
of  the  bay  was  covered  with  drifting  floes  of  young  ice  approxi¬ 
mately  one  inch  thick.  Strong  easterly  winds  with  gusts  to  50 
knots  completely  cleared  the  bay  of  ice  on  3  November.  Ice  began 
to  form  again  on  5  November,  and  a  ten-tenths  concentration  of 
young  ice  was  attained  by  7  November  with  exception  of  an  ice-free 
area  adjacent  to  the  pier. 

Commencing  8  October  and  terminating  7  November,  serial  tempera¬ 
ture  and  salinity  data  were  obtained  at  4  stations.  The  locations 
of  Stations  1  and  2  concurred  with  the  locations  of  Stations  1  and 
2  for  1959  as  shown  in  Figure  2.  Station  3  was  located  approximately 
100  feet  north  of  Station  1,  while  Station  4  was  located  just  off 
the  shoreward  end  of  the  pier.  The  data  arc  presented  in  Appendix 
VI,  Surface  temperature  and  salinity  values  for  Station  1  are  plot¬ 
ted  in  Figure  11. 

Data  were  taken  at  Station  4  for  comparison  of  the  oceanographic 
structure  outside  the  bubbled  area  with  that  of  the  water  column  at 
Station  1  during  the  early  period  of  ice  growth  on  the  bay.  Occupa¬ 
tion  of  Station  4  necessitated  breaking  through  the  ice  cover.  Siiarp 
rises  in  surface  water  temperature  were  observed  on  21  and  24  October 
easterly  winds  with  speed  maxlmums  of  51  and  48  knots,  respectively, 
were  recorded  on  these  dates.  Although  no  data  below  the  10-meter 
level  are  available,  it  is  evident,  as  indicated  by  the  temperature 
and  salinity  data  presented  in  Appendix  VI,  that  the  wind  affected 
vertical  mixing  throughout  North  Star  Bugt. 

On  10  October,  the  surface  temperature  at  Station  1  was  -1,54°  C 
the  salinity  was  32.30  °/oo.  On  15  October,  the  surface  temperatxire 
at  Station  2  was  -1.77°  C  with  grease  ice  forming  in  the  area;  sur¬ 
face  temperature  in  the  ice-free  bubbled  area  was  -1.68°  C.  Surface 
values  of  -1.8l°  C  and  32.82  °/oo  were  recorded  at  Station  1  on  29 
October;  the  bubbling  system  was  not  in  operation,  and  a  considerable 
amount  of  grease  and  slush  ice  was  forming  on  the  bay. 

Activation  of  the  bubbling  system  on  the  following  day  resulted 
in  quick  dispersal  of  all  ice  from  the  bubbled  area.  Light  grease 
and  slush  ice  being  swept  from  divergent  regions  and  transported 
beneath  the  surface  in  convergent  regions  confirmed  one  aspect  of 
the  proposed  model  of  induced  circulation.  Surface  temperature  of 
-1,82°  C  within  the  bubbled  area  indicates  supercooling,  since  the 
calculated  freezing  point  was  -1.79°  C. 


when  compared  to  data  obtained  at  Station  1,  those  obtained  at 
Station  4  on  5>  6,  and  7  November  indicate  that  vertical  trans¬ 
port  of  sensible  heat  was  not  a  factor  in  maintenance  of  the 
lee -free  area  adjacent  to  the  pier.  The  temperature  beneath 
the  lee  outside  the  bubbled  area  was  Identical  to  that  of  the 
Isothermal  water  column  at  Station  1. 

Subsequent  history  of  the  polynya  was  similar  to  that  of 
the  previous  winter. 
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TEMPCRATURErC) 


APPENDIX  III 


OCEANOGRAPHIC  DATA 
INNER  STATIONS,  1959 
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